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DYNAMIC CHARACTERISTICS OF A VIBRATING BEAM WITH 


PERIODIC VARIATION IN BENDING STIFFNESS 


JOHN S TOWNSEND 


ABSTRACT - A detailed dynamic analysis is performed of a 
vibrating beam with bending stiffness periodic in the spatial 
coordinate. Using a perturbation expansion technique the 
free vibration solution is obtained in a closed-form, and the 
effects of system parameters on beam response are explored. 

It is found that periodic stiffness acts to modulate the modal 
displacements from the characteristic shape of a simple sine 
wave. The results are verified by a finite element solution 
and through experimental testing. 
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DEFINITION OF SYMBOLS 


Symbol 


Def initon 


A - Peak stiffness variation, A = (Elmax - Elmin) / 2 

Anm - Factor defined by equation (A-2) 

dx - Differential element length 

d_ , ( ' ) - Denotes differentiation with respect to 

dx position 

Ela - Geometric average stiffness of beam design, 

Ela = (Elmax + Elmin) / 2 

EI(x) - Bending stiffness function 

k - Dummy variable referring to mode number 

L - Span length 

Lt - Period of the stiffness function 

m - Dummy variable referring to mode number 

rtf - Mass per unit length 

n - Vibration mode number 

Un(x) - Denotes nth mode eigenfunction, bending 

x - Distance along span measured relative to support 

Dimensionless horizontal coordinate, 'x = x/L 

Kronecker delta function, equation (10) 

Stiffness perturbation parameters, 6 = A / Ela 

Z - |4 2 


Ay 

X 




X 


n 


-A 

Dimensionless eigenvalue, \ = ^2 

Xn FT, 
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DEFINITION OF SYMBOLS (cont.) 


Symbol 


Def initon 


^ no 
i 

X nt. 
7T 

P 

o 

^ni 

0H2. 

Ct)r» 

J 

/ 

I 


Zeroth-order eigenvalue solution 
First-order eigenvalue solution 
Second-order eigenvalue solution 
- Pi 

Stiffness parameter, Q - 1,2,3 

ie., the number of half periods f) = 2 

Dimensionless Eigenfunction, = Un / L 

Zeroth-order eigenfunction solution 

First-order eigenfunction solution 

Second-order eigenfunction solution 

Denotes nth mode frequency, bending 

Partial derivative notation 

Integral sign 

Summation sign 


L / Lt 
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DYNAMIC CHARACTERISTICS OF A VIBRATING BEAM WITH 
PERIODIC VARIATION IN BENDING STIFFNESS 


INTRODUCTION 


Vibrating beam theory has been considered extensively 
in the literature for any number of variable property states, 
ranging from structures with changing cross-sectional geometry 
to those of a composite nature. Solutions are obtained either 
in closed-form for a few simple cases, or they are pursued 
using numerical techniques. In the present study, a perturba- 
tion expansion technique applicable to continuous systems is 
used to develop a closed-form solution to the problem of a 
vibrating beam with bending stiffness periodic in the spatial 
coordinate. Results are compared to a finite element solution 
and verified experimentally using forced vibration of a test 
span. To the knowledge of the author, this specific beam 
problem (static or dynamic) has not been addressed in the 
literature . 

Application of periodic stiffness is recognized in the field 
of vortex- induced motion of transmission power lines C13. 

In recent years a conductor, known as twisted-paired, has been 
developed that uses a variable diameter design to provide a 
changing conductor profile into the wind. Twisted-paired 
conductors are constructed by twisting together two identical 
standard round conductors with 360 degree twists occurring at 
set intervals along the span. The periodic nature of the twist 
causes a periodic variation in bending stiffness. Variable 
profile diameter results in non-uniform shedding of vortices, 
and hence excitation frequences, along the. span. Multiple 
vortex frequencies act to minimize wind energy transfer and 
detune vibration response. In conductor systems, influence 
of bending stiffness effects becomes extremely important in 
the vicinity of supports. 

The purpose of this report is to characterize the dynamic 
bending behavior of beams with periodic stiffness variation. 

Also, the models developed will provide insight into the behavior 
of similar type systems with changing property states. 
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EQUATIONS OF MOTION AND BASIC ASSUMPTIONS 


Consider the problem of the transverse vibrations of a 
straight beam with periodic variation in bending stiffness 
along its length. The beam is assumed to be simply supported 
and long compared to its cross-sectional dimensions, and dynamic 
shear distortions and rotary inertia are negligible. We will 
also make the usual simplifying assumptions that Hooke's Law 
holds and plane sections remain plane. Figure 1 shows a free 
body sketch of a differential element of the vibrating beam. 

We will proceed from the well-known differential equation 
of motion for the normal mode response of an undamped beam 



where EI(x) is the bending stiffness function, m is the mass per 
unit length, Un(x) is the nth normal mode displacement and CJ/j is 
the nth normal mode frequency. The bending stiffness function 
is given by 


EX (x) = £X^ + A cojx 22 

£-t 


where Ela is the geometric average stiffness of the beam 
design, Ela = CEImax + ElminI / 2 ; and A is the peak stiffness 
variation, A - CEImax - Elminl / 2. The period of the stiffness 
function is Lt. Figure 2 plots the function. Notice, at X = 0, 
the maximum flexural stiffness occurs, and the periodicity of the 
function is an even multiple of the span length. This particular 
stiffness function is characteristic of twisted-paired systems C13. 

An equation of motion that models an undamped, vibrating beam 
with periodic bending stiffness is determined by combining equations 
(1) and (2). The result is a fourth-order differential equation 
with variable coefficients. A closed-form approximate solution to 
this boundary value problem is obtained by using a variation of 
the Ray leigh-Schrodinger expansion C2,31. The solution in closed- 
form is extremely useful, since it clearly displays the influence 
of system parameters on response. Nayfeh C41 presents an applica- 
tion of a similar perturbation formulation for a simple linear 
second-order eigenvalue problem. 
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A dimensionless form of the governing motion equation is 



(See nomenclature section for a definition of terms.) 

and the corresponding boundary conditions for the case of simple 
supports are 


4n<-°) = “ 0 

and 

4>n (0) = <t>nO) = 0 


The quantity Cl + 6cosp7T xl is the dimensionless bending stiff- 
ness. For convenience, the tildes are dropped in the remaining 
analysis. The coefficient 6 is a measure of the magnitude of 
the stiffness variation, £ = A / Ela , and the parameter P is 
equal to the number of half periods of the stiffness function 
in a given span, P = 2 L / Lt. 


PERTURBATION EXPANSION SOLUTION 

The solution ( (f) , X ) of equation (3) is a function of the 
independent variable x and the parameters £ and P . If the 
parameter £ is equal to zero, the equation reduces to the case 
of a vibrating beam with uniform flexural stiffness whose eigen- 
functions and eigenvalues are given, respectively, by 


4'n ~ {I iu/n. nirx 


X 


Z 

n 




(7) 
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The coefficient is arbitrary. It is picked so that the eigen- 
functions are normalized according to the integral function. 



( 8 ) 


The above eigenf uctions are orthonormal; i.e.. 


I 

j ^ dx' - d m n 

C 

where, <5^-, the Kronecker delta function is specified as 


<5 


m n 


o m t n 

/ m = n 


( 10 ) 


When £ does not equal zero, equations (6) and (7) are no 
longer valid and corrections must be added to them. An approxi- 
mate solution is obtained by expanding both the eigenfunction and 
the square of the eigenvalue in the form of a power series in 
6 ; i.e. , 


0 


n 






(id 






6“ \n; 


( 12 ) 


where GD r . 0 and \ no are the eigenfunction and the square of the egen- 
value when £ equals zero; equations (6) and (7). An asymptotic 
expansion is generally valid only if £ is small. By definition, 
the parameter £ for beams with periodic stiffness may not be small, 
but it is always less than one. Corrections of the higher order 
terms are therfore negligible, and the series converges eventually 
to the correct solution. 
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Substituting equations (11) and (12) into equation (3) and 
equating like powers of 6 through order , we get the following 

system of differential equations: 

0 (€°) : <t>Z - Xn'.fino- 0 <131 


O(e') : 4>Z - Ko<Pn, - K$ no - £z [ <tL<&P*x ] < 34 . 


o(f) : C - X no 0 o ^X n ^„ o+ X n ,0 n| -^Mip7TX] <13. 


The problem of finding an approximate solution to equation (3) 
is now simplified to one of obtaining sequential solutions to 
equations (13), (14), and (15). To illustrate the procedure, the 
first-order correction is formulated in Appendix A. Equations (A-6) 
and ( A-8 ) define the correction terms of the eigenvalues and 
eigenfunctions . 

Using equation (A-8) and recalling that (p n = (p n0 + € (/) nt 
where P_=. 1,2,3, . . ., the general eigenfunction solution ^) n is 

expressed in terms of dimensionless parameters for the case when 
n / f) / 2 as 


= V 2 lUl tiitx 


€ 77 

V? 


€ 71 



~{2 (n-p) A -n*\ 


( n +p)?rx 




If n = p/ 2 , the last term in the series is secular and hence 
vanishes . 
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Using equations (7), (12), and (A-6), the eigenvalue 
solution is given in terms of the first-order correction. 
The result is expressed in terms of the dimensionless variables 
as 

X n = r? 2 7T 2 (i - -j Y <17 


Equation (17) is valid only for the vibration mode n = yO / 2 . 

For all other vibration states the first-order correction term 
is zero. The second-order perturbation solution of the eigenvalue 
is determined using the same techniques previously developed for 
finding the first-order terms. Details of the formulation are out- 
lined in Reference C1D. The general eigenvalue solution \ of 
the second-order expansion is given as /xr ' 


X 


n 


1 2 . 

n n 


/ _ L -<t 

(n+p) A 1 

€ 

1 

5 

i 

cr 
\ 


1 2 4 

_(n+p)*-n\ 

4 

(-n-pf-rf _ 



The first-order correction term in equation (18) is equal to zero 
for the vibration modes where n /p/ 2. For the case defined by 
n = P/ 2, the last term in the second-order correction is specified 
to vanish (ie., this term is secular from ^P'n I solution). Notice 
that the eigenvalue simplifies to the case of a beam with 

uniform stiffness when the perturbation parameter 6 is equal to 
zero. 


The general behavior of a beam with periodic bending stiff- 
ness variation is given by equations (16) and (18). As the 
vibration state approaches the anomaly occurring at n = yO / 2, the 
eigenfunction solution deviates from a simple sine wave displace- 
ment curve to a mode shape comprising other harmonics. For the 
case when n = £)/2, the eigenfunction returns to the sine wave shape. 
The eigenvalue solution responds in a similar nature. At n =C / 2 
a jump in the eigenvalue occurs, since for this mode the harmonics 
of the stiffness function are secular. Results of the closed-form 
perturbation analysis have been checked using finite element results 
and results of experimentation. 
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EIGENFUNCTIONS AND EIGENVALUES 


Effects of the perturbation terms on the eigenfunction 
solution are exemplified in Figure 3. The magnitudes of the 
bracketed terms in equation (16) are plotted as a function of 
the p/n ratio. Two distinct ranges are apparent; P In < 1 
and near D / n = 2. In these ranges, the perturbation effects 
are the strongest. At pi n = 2, a vibration state is defined 
where the lengths of the vibration loops match the period of 
the stiffness function. Based on the stiffness definition given, 
this mode defines maximum stiffness at the nodes of the vibra- 
tion loops and minimum siffness at the antinodes. At pin = 1, 
both the maximum and minimum stiffnesses occur at the node 
positions in alternating sequence along the span. 

Figures 4 and 5 give eigenfunction comparisons of vibration 
modes in the general modal solution range. The vibration displace- 
ment amplitudes are normalized and plotted verses the normalized 
horizontal span coordinate (X/L) . Recall, boundary conditions are 
simple support and beam orientation at the supports is for maximum 
stiffness. Figure 4 plots eigenfunction solutions near the anomaly 
pin = 2, where P = 64 and € = 0.4. The effect of periodic 
stiffness is to modulate the displacements of those vibration modes 
approaching the anomaly at n = /0 / 2, or for this case mode 32. 

Similar displacement curves as patterned for modes 31 and 33 are 
characteristic of all vibration modes near the anomaly. Close exam- 
ination of the eigenfunctions reveals that the node (or antinode) 
locations are adjusting themselves along the span, and the longer 
vibrating loops result in lower midloop displacement amplitudes. 
Apparently, the beam attempts to minimize the elastic strain energy 
stored within the dynamic span by adjusting the lengths of the 
vibrating loops until the same average bending siffness exists across 
each individual loop. Equalizing the loop stiffnesses may require 
the loops to have different legnths depending on the vibration mode, 
and a longer loop has greater mass. An equal partioning of potential 
energy and thus kinetic energy between each of the loops results in 
lower vibration amplitudes for the longer vibrating loops. Loop 
stiffness calculations verify this reasoning. 

At the anomaly, modulation in the mode shape disappears, since 
for this case the individual loop stiffnesses are equal (ie., the 
lengths of the vibration loops match the period of the stiffness 
function). The same basic reasoning holds true for the case where 
P /n = 1, see Figure 5. Here also, the mode shape is sinusoidal 
- no modulation; and average loop stiffnesses are equivalent with 
maximum and minimum values defining the nodes of each loop. Add- 
itional cases identified in Figure 5 are for small and large values 
of p/n. As pin approaches zero, perturbations in the mode shape 
increase. At the opposite end of the scale, where the ( Q In para- 
meter goes to infinity, the displacement shape is sinusoidal. 

Figure 6 characterizes the eigenfunction solution as a function of 
the perturbation parameter, £ . As the magnitude of Epsilon (ie., 
the stiffness variation) increases, the modulation effects become 
more pronounced. 
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FIGURE 3 EFFECTS OF PERTURBATION TERMS ON EIGENFUNCTION 
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IGENFUNCTION SOLUTION 



FIGURE 5 EIGENFUNCTION SOLUTIONS IN CHARACTERISTIC RANGES OF VIBRATION 









OF P 




Eigenvalue solutions of the 2nd order perturbation expansion 
for different cases of P are plotted in Figure 7 versus the respec- 
tive vibration mode. At the anomaly, f) I n = 2, a jump in the eigen- 
value occurs. The intensity of the jump increases with p . Physic- 
ally, the jump identifies a sharp change in frequency (or stiffness) 
between characteristic modes of vibration. Mathmatically , the jump 
is equivalent to removing the secular nature of the stiffness func- 
tion from the eigenfunction solution. 

The closed-form perturbation expansion solution has been 
verified through comparisons with a finite element solution C1D. 
Although the findings are not formally documented herein, 
agreement between the analytical results is excellent. Some 
discrepancy does occur in the vicinity of the stiffness anomaly. 

This is apparently due to a sudden change in bending stiffness. 
Nevertheless, the qualitative picture remains the same. Two 
characteristic effects of periodic bending stiffness on dynamic 
response are determined: (1) periodic stiffness forces an anomaly 

in the system which results in a jump in the natural frequency, 
and (2) periodic stiffness acts to modulate the modal displacements 
in distinctive ranges of f)/ n. A qualitative explanation of the 
modulation and its effects on beam response is given in terms of 
energy principles. 


EXPERIMENTAL VERIFICATION 


A series of tests were designed to investigate the dynamic 
response of beam type systems which have a periodic variation in 
bending stiffness. A stiff-string structure, known as a twisted- 
paired conductor, was the test candidate in the program. 

Table 1 summarizes the experimental test parameters and Figure 8 
shows a photograph of the test span. Periodic variation in 
diameter profile of the twisted-paired conductor is compared to 
the uniform diameter of a standard conductor design. In stiff - 
string systems, elastic strain energy is stored in tension and 
bending. If tension is constant along the span, then tension has 
minimal effect towards equalizing the variable flexural stiffness 
of the vibrating loops C1D. In other words, tension effects do 
not mask the effects of stiffness variation. 

Experimental data are compared to the finite element results 
for free vibration since the fixed boundary conditions are applic- 
able. This type of boundary support keeps the end losses to a 
minimum. Internal damping of the conductor was also reduced by 
applying a high tension line force. The procedure of minimizing 
conductor system damping is necessary; higher harmonics are 
difficult to excite if mechanical damping is significant. The 
testing program used forced vibration response to study free 
vibration. If the span is tuned properly to a single natural 
frequency, contributions from all other harmonics are minimal. 

The vibration exciter unit was positioned near the span center 
to eliminate even harmonics from the general response. The added 
mass of the moving shaker element and span attachment fixture 
resulted in a shortening of the drive loop, and thus a lowering 
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FIGURE 8 PHOTOGRAPH OF TEST SPAN, TWISTED-PAIRED CONDUCTOR 




MODE 



FIGURE 10 COMPARISON OF EIGENSOLUTION WITH 
EXPERIMENTAL TEST VALVES, MODE 19 




of its vibration amplitude. No attempt was made to decouple 
shaker mass from the conductor span. A V- scope attached to the 
center of each vibrating loop was used to measure the midloop 
amplitude displacements. The device is inexpensive and its accuracy 
is remarkable at 0.01 inches. 

Typical test results are presented in Figures 9 and 10. 

Finite element displacement amplitudes are normalized using a 
method previously outlined, see equation 8. Experimental 
amplitudes are normalized to one of the measured values - chosen 
in arbitrary fashion. Span length is used to nondimensionalize 
the horizontal coordinate. Since mode shape is symmetric about 
the span center, data results are shown only for half the span. 
Figure 9 gives the comparison for mode 27, an eigenfunction near 
the anomaly in the system occurring at n = 32. Although the com- 
parison is not exact, the modulation in the eigenfunction response 
is proved physically to exist. The same general results are re- 
ported for all other modes near the anomaly. Figure 10 compares 
the experimental and analytical data of mode 19, an eigenfunction 
well removed from the stiffness anomaly. Agreement between the 
analytical and test data is excellent and the mode shape is 
sinusoidal. Some discrepancy does occur in the node positions 
near the drive location where the measured loop lengths are shorter 
Shacker attachment changes the stiffness and mass of the drive loop 
Shifting of the nodes tends to compensate for these effects. 


CONCLUSIONS 


The dynamic response of beams with periodic stiffness 
contrasts significantly with the vibration behavior of standard 
beams. Linear vibration theory was used to develop a stiffness 
model and characterize response. Using a perturbation expansion, 
a closed-form solution of free vibration was formulated for the 
case of simple supports and periodic stiffness variation. The 
technique worked exceptionally well when the stiffness parameter 
was slowly varying. Applications of variable tension, mass, and 
area are natural extensions of the theory. The main conclusions 
are summarized below. 

1. Periodic bending stiffness forces an anomaly in 
the system which corresponds to the vibration state 
where the loop length matches the period of the 
stiffness function. Physically, the anomaly denotes 
the vibration mode for which loop stiffness changes 
most rapidly. The result is a jump in natural freq- 
uency. The perturbation solution loses some accuracy 
for those vibration modes near the anomaly; however, 
the qualitative characteristics of the response 
remain the same . 

2. The stiffness parameter acts to modulate the modal 
displacements in two distinct ranges of vibration: 

p /n < 1 and near f) / n = 2. Experimental evidence 
is presented which supports these findings. 
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3. Modulation in modal displacements is explained 

in terms of energy principles. The beam attempts 
to minimize the elastic strain energy stored 
within a dynamic span by adjusting the lengths of 
the vibrating loops until the same average bending 
stiffness exists across each individual loop. 
Equalizing the loop stiffnesses may require the 
loops to have different lengths depending on the 
vibration mode, and a longer loop has greater mass. 
An equal partitioning of potential energy and thus 
kinetic energy between each of the loops results in 
lower vibration amplitudes for the longer vibrating 
loops. 
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APPENDIX A 


FIRST-ORDER CORRECTION TO THE PERTURBATION SOLUTION 


The first step in the first-order correction to the perturbation 
solution is to substitute the zeroth-order solution, equations (6) 
and (7), into equation (14). After taking the appropriate derivatives 
and using trigometric identities, the result simplifies to 

^ A A [ 

<£„, - n 7T <£ n , = y 2 ( X n , iwv nnx 


( A-l ) 

- [ ( n +P) 7rx] - ^{r)-p)tijrJ[{r)-P)TTX] • 

4 

Next, assume that the solution 0 ni can be expressed as a linear com- 
bination of the zeroth-order eigenfunctions . ; 


0 n , = Z Anm <**?> 

i 

This solution satifsies the boundary conditions on : fir\\ • Taking 
derivatives and substituting into equation (A-l) yields 

£ 7T 1 (m 4 -)?' , )Anm NfVj m7TX = 

Xni Avv'-nTTA - iLir (n+pfw^+P^X <a-3> 

2- ^ 2. 

- (r?-p) -P) / ttX 

Multiplying equation (A-3) by sin k/Tx and interating from 0 to 1 
using the orthonormal property, equation (9), we obtain 
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r 

7 T* [ K - n* ] AnK = Xm <$m -ftSr *(n + pf J M^{n+p)lTX^K 7 Txd){ 


I 

— v ^{n-p) JHu^( n-p)TTX Ma^-> KFX d X 


(A-4) 


If k = n r the left-hand side of equation (A-4) vanishes, hence 

I 

X rtl = n\ A (n+pf J Awv,(r)+p) 7 rX^i 77 TX<ix 

o 

l 

+ r?V(>?-pf J k^(ri -p)nX!^v'TTxdK 


( A- 5 ) 


The above integral expressions evaluate only when n = v/2, that 

is when the vibration mode number corresponds to the span length- 
stiffness function ratio (n = L/Lt). This is the anomaly that 
makes the periodic stiffness problem so interesting. 


Equation (A-5) then calculates the eigenvalues of the first- 
order expansion as 


X 


hi 


- J?*7T 4 


U>HeN T) = p I 2 


(A-6) 


Note, Xm = 0 for all other values of n. The above condition removes 
the secular terms from the solution when k = n. If k ^ n, then 
equation (A-4) simplifies to 


A nK ~ 


I 

2 - C 

J hU^(n+P)TTK KttXcU 


( A-7 ) 


I 


^-LR P) f r) - p)ttx mavj KnXclx 


z 



Because of the stipulations on the parameter P , Ank in equation 
( A-7 ) calculates non-zero values only for the two cases: 

( 1 ) k » n + p and ( 2 ) k = n - p , wher n / (0 /2. The 

general solution of the eigenfunction <p A( , given by equation (A-2), 

is finally expressed as 

4 rT r l+P)'JTX 

VI l(n+P) -n*\ 

(A-8) 

{z [(v-py- n*] ^ j(n " P ) 7rX + JIA m b^rmX 

For the vibration mode corresponding to n = p/2, the second term 
in equation (A-8) vanishes, since for this mode it is secular in 
nature. Keeping this in mind we can say that the m )/ solution is 
valid for all vibration states where P= 1,2,3,. . . The 
coefficient Ann is determined by the normilizing function 


/ 

^ ^ dX ~ 0 

a 


(A-9) 


For this case, Ann is calculated to equal zero. 
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